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In the crystal of K 3 H(S0 4 ) 2 , or K 3 D(S0 4 ) 2 , dimers (S0 4 ) • • -H- • -(S0 4 ), or 
(S0 4 ) • • -D- ■ -(S0 4 ), are linked by strong centrosymmetric hydrogen or deuterium 
43 bonds whose O- • -O length is ~ 2.50 A. We address two open questions, (i) Are H 



or D sites split or not? (ii) Is there any structural counterpart to the phase tran- 
sition observed for K 3 D(S0 4 ) 2 at T c « 85.5 K, that does not exist for K 3 H(S0 4 ) 2 . 
Neutron diffraction by single-crystals at cryogenic or room temperature reveals no 
structural transition and no resolvable splitting of H or D sites. However, the width of 

O ■ 

the probability densities suggest unresolved splitting of the wavefunctions suggesting 
rigid entities Hxi/ 2 — Hri/ 2 or Dxi/ 2 — D#i/ 2 whose separation lengths are /h ~ 0.16 A 
or Id w 0.25 A. The vibrational eigenstates for the center of mass of Hii/ 2 — H^ 4 / 2 
revealed by inelastic neutron scattering are amenable to a square-well and we sup- 



0> pose the same potential holds for Dli/ 2 — D#i/ 2 . In order to explain dielectric and 

calorimetric measurements of mixed crystals K 3 D(!_p)H p (S0 4 ) 2 (0 < p < 1), we re- 
place the classical notion of order-disorder by the quantum notion of discernible (eg 
Dli/2— D/ji/2) or indiscernible (eg H £1 / 2 — H R1 / 2 ) components depending on the sep- 

h: 

aration length of the split wavefunction. The discernible-indiscernible isostructural 
transition at finite temperatures is induced by a thermal pure quantum state or at 
K by p. 

PACS numbers: 61.66.-f, 61.05.F-, 64.70.Tg 
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I. INTRODUCTION 



Symmetric hydrogen bonds O- ■ -H- • -O are of fundamental importance to many disciplines 
across physics, chemistry, or biology, in the gaseous, liquid or solid state. Their properties 
depend on essentially the O- • -O bond length (say Roo), such that a shorter -Roo corre- 
sponds to a stronger bonding and a lower frequency for the OH stretching mode (VOH). 
This dependence on i?oo accounts for temperature, pressure, and isotope effects. As Roo 
decreases, the potential operator for the stretching coordinate is expected to evolve from a 
double-well with a high barrier (compared to kT) to a single-well for the shortest bonds. The 
latter is widely thought of as mimicking symmetric intermediates for proton transfer along 
the vOH coordinate in (bio) chemical reactions. This is an important field of investigations 
for computational chemistry and our purpose is to determine effective potentials based upon 
experiments for benchmark tests of theoretical methods. 

In crystals, strong symmetric hydrogen bonds are encountered for small entities, such as 
chelates or dimers, well separated from each other. Bond lengths, symmetry and proba- 
bility densities at proton sites are best characterized with single-crystal neutron diffraction 
(SXND), compared to x-ray (SXXD). Inelastic neutron scattering (INS) reveals OH tran- 
sitions with spectacular contrast of intensity, compared to infrared or Raman. Probabil- 
ity densities in the momentum space can be determined with neutron Compton scattering 
(NCS). 

The best characterized short symmetric hydrogen bond in a crystal is that of the chelate of 
the potassiumhydrogenmaleate, KH(OOC-CH=CH-COO), for which R 00 = 2.427(1) A at 
20 K.- This is one of the shortest distance ever reported. INS spectra reveal numerous z/OH 
transitions consistent with a funnel-shaped potential, such that the bare proton is localized at 
the center in the ground state, or largely delocalized over the oxygen orbitals in the excited 
states beyond ~ 500 cm -1 . Similar potentials have been proposed for centrosymmetric 
dimers of potassiumhydrogenbistrifluoroacetate, KH(CF3COO)2 and the cesium analogue, 
CsH(CF 3 COO) 2 , both with .Roo = 2.436(4) A at cryogenic temperatures.- In crystals, the 
existence of symmetric double-wells for longer bonds is controversial. For the crystal of 
potassiumdihydrogenphosphate (KDP, Roo ~ 2.50 A), Reiter et al.- claimed that NCS 
suggests a symmetric double-well in the paraelectric phase above T c rj 124 K, or a single 
minimum in the ferroelectric phase below T c . However, this is clearly excluded by the crystal 
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symmetry and vibrational spectra.- Quasi-symmetric double wells have been reported for 
centrosymmetric dimers of the hydrogencarbonate family (MHCO3 with M = K, Rb, Cs),-£ 
or benzoic acid,- 1 ^ for which Roo ~ 2.60 A, but these hydrogen bonds are not symmetric. 
Apart from crystals, the chelate of the malonaldehyde (3-hydroxy-2-propenal) molecule in 
the gas phase (-R00 ~ 2.57 A) is the best characterized symmetric double- well.- 1 ^ 

In the present work we report experimental and theoretical studies of the crystals of 
trispotassiumhydrogendisulfate, K 3 H(S04)2, and of the isostructural K 3 D(SC>4)2.— ~— The 
unit cell of K 3 H(S0 4 ) 2 is comprised of 4 equivalent centrosymmetric dimers, S0 4 • • -H- • -S0 4 
for which #00 « 2.515(2) A at 293 K or « 2.487(2) A at 100 K. These small dimers are well 
separate by K nuclei. The monoclinic group symmetry is A2/a and there is no structural 
phase-transition detectable by SXXD from 4 K to 300 K. Noda et aL 12 suggested that split 
proton sites at room temperature turn into single proton sites at low temperature, as Roo 
becomes shorter (geometric effect). INS reveals OH bending modes at 1250 or 1600 cm -1 
and the z/OH mode at rj 57 cm -1 . This is the lowest frequency ever reported for any 
hydrogen bond in a crystal. 14 This could be an evidence to suggest an asymmetric double- 
well with a low barrier, but this is forbidden by the crystal symmetry. We report below 
SXND measurements showing that a double-well is unlikely at every temperature and we 
propose a new potential operator for INS. 

Several members of the M 3 H(X04)2 family (M = K, Rb, Cs; X = S, Se) demonstrate 
remarkable isotope effects which reveal intriguing differences between H and D bonds. In the 
specific case of K 3 D(SC>4)2, calorimetric and dielectric measurements reveal anomalies at- 
tributed to a phase transition at T c 85.5 K, whereas there is no transition for K 3 H(S04)2-— 
So far, T c is conceived of as an order-disorder transition from an antiferroelectric phase 
(A2) below T c to a paraelectric phase (A2/a) above T c ,— ~— A widespread explanation is 
that deuterons are delocalized in a symmetric double-well above T c and statistically local- 
ized below T c .— ii£>il"i2 Two models have been proposed to explain the isotope effect, (i) 
The pseudospin Ising model with markedly different tunnel splitting for H or D.— liLiS (ii) 
Markedly different geometric effects for H or D bonds.— iI2J&22 Both models are weakened 
by the lack of unquestionable evidence to suggest A2 symmetry consistent with asymmetric 
O- ■ -D- ■ -O bonds below T c . We have, therefore, carried out SXND measurements in order 
to observe whether or not there is a structural counterpart (possibly A2 < — > A2/a) to the 
dielectric and calorimetric anomalies. 
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This article is organized as follows. In Sec. [TTJ we present SXND data for K 3 H(S04)2 and 
K 3 D(S04)2 at cryogenic and room temperatures. It transpires that the space group is A2/a 
at every temperature for both crystals and the best refined structures are obtained with H 
or D at the center of symmetry. Statistical distributions and antiferroelectricity below T c 
are, therefore, ruled out. In Sec. 1111} we show that SXND and INS are together consistent 
with unresolved split protons in a square-well and the same model is supposed to hold for 
deuterons. In Sec. IIVI we propose a new interpretation within the framework of quantum 
mechanics of the dielectric and calorimetric anomalies for KD( 1 _ p )II o (S04) 2 , < p < 1. We 
emphasize that these crystals are macroscopically quantum, to the least, below T c . In Sec. 
|V]we show that the same model holds for the isostructural Rb 3 H(S04) 2 . 

II. NEUTRON DIFFRACTION 

Single- crystals were obtained by slow cooling of H 2 S04(D 2 S04);1.2 K 2 S0 4 in H 2 0(D 2 0). 
Prismatic colorless specimens wrapped in aluminum were loaded in a closed-cycle-refrigerator 
whose temperature was controlled to ±1 K. Data were collected with the four-circle diffrac- 
tometer 5C2 at the Orphee reactor (Laboratoire Leon Brillouin).— Every structural param- 
eter computed with CRYSTALS 22 was allowed to vary independently.— 

Inspection of intensities for absent reflections confirms the monoclinic space group A2 / a 
for both K 3 H(S04) 2 and K 3 D(S04) 2 at every temperature (see Tabled]). There is no visible 
symmetry breaking to suggest a statistical distribution of H or D. Best refinements were 
obtained with H or D at special positions (Tables [Til and IHIj) and every attempt to split 
these sites was less satisfactory. The best refined structures are comprised of Z = 4 dimer 
entities per unit cell. The slightly distorted S0 4 tetrahedra are consistent with spontaneous 
dielectric polarization. They are linked through centrosymmetric O- • -H- • -O or O- • -D- • -O 
bonds practically parallel to (a, b) and at ~ ±30° with respect to b (Fig. [TJ). The supposedly 
antiferroelectric phase of K 3 D(S04) 2 below T c is not confirmed and both the Ising model 
and geometric effects are rejected. 

Unit cell parameters (Table HJ, positional parameters (Tables ILT1 and ITTTj) ). and thermal pa- 
rameters for heavy nuclei (Tables lIVl andlVl) are similar for both isotopomers but temperature 
effects are different. For K 3 H(S0 4 ) 2 , the volume of the unit cell increases by ~ 3% between 5 
K and 293 K, R ouo increases from 2.463(1) A to 2.496(1) A (Ai? OH o = 0.033(2) A), whereas 
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K 3 H(S0 4 ) 2 



K 3 D(S0 4 ) 2 



FIG. 1. Projections of dimer entities onto (a, b) planes. Thermal ellipsoids represent 50% of the 
probability density for nuclei. Sulfur is yellow, oxygen is red, potassium is black, hydrogen or 
deuterium are grey. 

the S— O bond lengths (i2|o) are practically unchanged (Table IVT1) . For K 3 D(S0 4 ) 2 , the unit 
cell is practically temperature independent, -Rodo shortens very slightly from 2.532(1) A to 
2.520(2) A (AifoDO = -0.012(3) A) and the R^s increase significantly, what is in accor- 
dance with a change of the spontaneous dielectric polarization. 

The thermal parameters for H and D are markedly anisotropic. The diagonal factors 
U X (T) and U Z (T) for coordinates x or z in (a, b) planes, parallel or perpendicular to O- • -O, 
respectively, are given in Table IVTT1 The corresponding half- widths at half- height (HWHH), 
A Q (T) (a = x, y, z), of the gaussian profiles 



P(a,T) 



exp[-a 2 /U a (T)} 



(1) 



are given in Table IVIIII Because the rather modest thermal parameters for heavy nuclei 
are virtually identical for both isotopomers, U a (T) and A a (T) are largely representative 
of vibrational displacements. For K 3 H(S0 4 ) 2 , U^(5) and £7^(293) suggest a markedly an- 
harmonic z/OH mode at a low frequency and/or an unresolved splitting. For K 3 D(S0 4 )2, 
[/J? (14) S> U^(12) is at variance with the mass effect anticipated for usual single or double- 
well potentials. 
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FIG. 2. INS spectra of K 3 H(S0 4 ) 2 and Rb 3 H(S0 4 ) 2 crystal powders at 20 K, after ref. [l4]. 
bands at « 600 cm -1 were assigned to SO4 entities. 

III. THE INFINITE SQUARE- WELL 



The 



Since the double-well proposed in Ref. [14J is rejected by SXND, we seek a single- 
well potential for a proton oscillator (m = 1 amu) consistent with the z/OH transition 
observed at (57 ± 3) cm -1 and with the probability density profile. A (quasi)harmonic 
model is inappropriate because the mean-square amplitude in the ground state, (x 2 ) ~ 0.29 
A 2 . 21 should be one order of magnitude greater than Uf(5). Consequently, we envisage 
another textbook potential that is the infinitely deep square-well such that V = for 
—an/2 < x < an/2 and V = 00 otherwise.— The eigenstates are: 

n 2 h 2 

(2) 



ima 




11 



2 rniix + an/2) 

sin 

a H a H 



n > 1. 



The zero-point energy is Ef = (Ef - E?)/3 = (19 ± 1) cm" 1 and a H = (0.47 ± 0.01) A. 
Table [IX] shows that E* - Ef = {n 2 - l)Ef, 2 < n < 6, compare favorably with INS 
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transitions (Fig. [2]). This is an encouragement to pursue this model. The eigenfunctions 
are represented in Fig. [3] and the expected probability density, 




n — 1 — i — 1 — i — 1 — i — 1 — r 

-0.2 -0.1 0.0 0.1 0.2 



x(A) 

FIG. 3. Eigenstates of the square- well along the OH stretching coordinate 



ty 2 (T) = ^ l ^ v v — 17/ /|yt| (3) 

1 ) E l ^M-(E l -E 1 )/kT) ' 1 ) 

can be confronted with (pQ). At 5 K, the leading term \ipi\ 2 gives a bell-shaped profile whose 
HWHH is Af x = a n /A « 0.12 A. This is significantly less than A*(5) « 0.156 A. We 
attribute the difference to a splitting of the wavefunction, what means that the energy levels 
refer to the center of mass of a rigid entity, say Hli/ 2 — Hm/ 2 , whose separation length is Zh 
and whose probability density is the sum of two unresolved profiles \ipi{x ± /h/2)| 2 /2. Our 
best estimate is Zh = (0.16 ± 0.01) A (see Fig. Hj). Then, \l/ 2 (293) computed with the same 
parameters cannot be compared directly to a gaussian profile. The best we can conclude is 
that the HWHH corresponding to the plateau at half-height, namely (0.20 ± 0.05) A, is in 
qualitative agreement with A^(293). At this modest level of precision, there is no visible 
change of an or /h consistent with the expansion of the O- • -O bond. There is no visible 
temperature induced geometric effect. 

For K 3 D(S04)2, OD transitions are unknown. We tentatively suppose ao ~ an, ~ 
E*/2, and |^°| 2 « | 2 - Then, A°(14) is consistent with Z D w 0.25 A (Fig. EJ. (Note 
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T = 293 K 




A 

FIG. 4. Solid lines: Computed probability densities along the OH stretching coordinate at 5 K 
(bottom) and 293 K (top) for a rigid H4/2 — H4/2 entity. The length is Zh = 0.16 A and the center 
of mass experiences a square-well whose width is 0.47 A. Dashed line: probability density for each 
half proton. +: Gaussian profile whose variance is Uf(T)/2 (see Table EDJ. 



that because of the mass effect, \l/ 2 (14) includes a minor contribution of \ip2{x ± fo/2)| 2 /2 
to ~ 5% of the total density). At 293 K, the HWHH of the calculated profile is roughly on 
the order of A°(293). Note that Z D - l B « 0.09 A is rather close to i? DO - Rouo ~ 0.07 A 
at low temperatures. This suggests a geometric effect upon isotope substitution. 

So far, the analysis of the density profiles is reasonably well established for K 3 H(S0 4 ) 2 at 
low temperature, since and Zr are determined by different measurements. For K3D (804)2 
at low temperature, we ignore whether the assumption fh is correct, so Z D should 

be regarded with caution. In addition, it is impossible to conclude whether or not the model 
is pertinent at elevated temperatures. These drawbacks are wiped away in the next section 
where we show that dielectric and calorimetric data reported for mixed isotope crystals 
K 3 D( 1 _ p )H p (S04)2 are precisely represented by this model. 
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A 

FIG. 5. Solid lines: Computed probability densities along the OD stretching coordinate at 14 K, 
85 K, and 293 K for a rigid D L1 / 2 — Dki/2 entity. The length is Id = 0.25 A and the center of mass 
experiences a box whose width is 0.47 A. Dashed line: contributions of the n = 1 and n = 2 states 
to the probability density. +: Gaussian profile whose variance is U^(T)/2 (Table IVTH) . 

IV. ISOTOPE EFFECTS AND PHASE TRANSITIONS 

Because centrosymmetric dimers have no electric dipole, there is no significant coupling 
terms between proton or deuteron vibrations which are, therefore, degenerate. This is 
confirmed by INS. The Ising model is, therefore, irrelevant. 

Consider an ideal defect-free crystal of K 3 D(S0 4 ) 2 represented by a periodic lattice in 
3-dimension of N indistinguishable rigid entities, T>Li/2j— Dm/2.?, centered at every node j. 
M is on the order of Avogadro's constant. The centrosymmetric wavefunction for the center 
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of mass Xj is 



lf> nj (Xj, l D ) = -^=[lp n (Xj + In/ 2) - ipnixj - Z D /2)], (4) 



The lattice can be represented by TV Bloch-wavefunctions indexed r (1 < r < A/"), 

1 N 

^ n (:r,Z D ,k nr ) = jj : ^2tp nj (x j ,l D )expik nr .L j , (5) 

i=i 

where x is a nonlocal observable independent of j, k nr is a wave vector and Lj is a lattice 
vector. Bloch-states are degenerate with respect to k nr , so the interpretation presented 
below is k nr independent. For K 3 H(S04)2, Bloch-states should be antisymmetrized if protons 
effectively behave as fermion, but this would be of no consequence to our interpretation. 

It is worth emphasizing that a symmetric double- well giving rise to tunneling is excluded 
because there is no symmetry plane perpendicular to O- ■ -O compatible with a symmetric 
ground state 2~ 1 / 2 [|L) + |-R)].— Likewise, symmetric displacements of £>Li/2j with respect 
to D^!/ 2 j are forbidden. Only displacements of the center of mass of the rigid entity are 
allowed. 

Below T c , K 3 H(S04)2 and K 3 D(S04)2 demonstrate different probability densities for H 
or D, different S— O bond lengths, opposite temperature effects for their dielectric constants 
(de H /dT < while de^/dT > 0) and e D reaches a maximum at T = T c . Above T c , probabil- 
ity densities, bond lengths and temperature effects for the dielectric constants are similar: 
de R /dT < 0; de° jdT < 0. From these observations we tentatively infer that the S— O bond 
lengths, hence the dielectric constants, could correlate with the splitting of the wavefunc- 
tion or, more specifically, with the nonlocal "indiscernibility degree" in the ground state 
proportional to the overlap defined on a — 1 scale as: 



(o-Q/2 

2 f . tt(x + a/2 + 1/2) . vr(x + a/2 - 1/2) , 

ZiUJ = - / sin sin dx 

a J a a 

-(0-O/2 
a — I nl 1 a — I 

= cos 1 — sin7r 

a an a 



(6) 



Then, we tentatively set the discernibility-indiscernibility boundary (DIB) to Ii(ldi) = 1/2, 
for Idi = 0.183 A. Consequently, Xi(Zh) > 1/2 at every temperature (see Fig. [6] a) and 
the S— O bonds are virtually temperature independent. In contrast, Zi(fo) < 1/2 at 14 K 
and we assume that indiscernibility emerges at elevated temperatures from a thermal pure 
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quantum state leading to coherent oscillations of the center of mass: 



|* D (t)) 



Zo%\il> n )exp(iE%t/h) 



n>l 



a 



D|2 



ex V [-(E»-E?)/kT} 



> . 



(7) 



Eexp[-(^ D -^i D )AT] 

n>l 

The temperature law for the indiscernibility degree is then 

X D (T) = |a?i 2 ^ D ) + J2 E ^ m \a»\ 2 , (8) 

n>l m>n 

where p n>rn is the indiscernibility degree of the superposition \ri) + e l<p \m) for which the mean 
position of the center of mass is:— 



COS 



An 



x sm sin ax 



[(n 



m 7r 



[cos(n — m)7r — 1] 



[(n + m)7r] 2 



[cos(n + m)iT — 1] 



(9) 



If m 7^ n + 1, X n m sa and ^ n>m ~ 0. Otherwise, X n>n+ i rs — 2a7r 2 and the mean amplitude 
of oscillations, (/) = 2y / 2a7r~ 2 , determines the time- averaged mean-separation of the split- 
wave: (Z D ) = l D - (I) = (0.11 ± 0.01) A. Then, (3 n<n+1 = Zi((Z D » = 0.77 ± 0.03. The DIB 
X D (T di ) = 1/2 gives T di = (86 ±10) K, in accordance with T c . Our assumptions Ti{l di ) = 1/2 
and E® ~ E^ /2, our estimate Zd = 0.25 A and the thermal pure quantum state (jTj), are thus 
validated at the macroscopic level of the phase transition. We conclude that the square-well 
is largely isotope and temperature independent, to the least below T di . 

Our interpretation of the isotope effect is as follows. For K 3 H(S04)2, Ii(Zd) > 1/2 
at every temperature, the S— O bonds are practically temperature independent, and the 
monotonous decrease of e H upon rising T is due to thermal fluctuations of the crystal lattice. 
For K 3 D(S0 4 ) 2 , Zi(fo) < 1/2 for T < T di , the S-0 bonds are longer and e D < e H . As T 
increases, Z D (T) increases, the S— O bonds shorten, e D increases and reaches a maximum at 
= T c . Since the probability density computed at 85 K is similar to that at 293 K (see 
Fig. EJ, we infer that the bond lengths are practically temperature independent above Tdi, 
so the monotonous decrease of e D is analogous to that of e H . 

So far, the cusp at T = T di for e D reveals a discontinuity that is not accounted for by dSJ). 
Consequently, we suppose that the phase transition is actually triggered by a discontinuity 
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"/(A) /(A) 

FIG. 6. (a) The indiscernibility degree of the ground state, T\(l), as a function of the separation 
of the split-wave. Z H = 0.16 A, X x = 0.60 for K 3 H(S0 4 ) 2 - h = 0.25 A, X 1 = 0.27 for K 3 D(S0 4 )2. 
Idi = 0.183 A, Xxdi = 0.50 at the discernible-indiscernible boundary, (b) The transition parameter 
TV =Xiill< l di , TV = 1 if I > l di . 



of the indiscernibility parameter, say XV\, such that: (i) XV i — X\ if I > Z^i or = 1 if 
/ < J* (see Fig. Mb); (ii) ZP D (Z) = Z D (T) and /3 n , n+1 = Zi«fo)), if Z < Z di , or ZP D (Z) = 1 
and ^ n ,n+i = 1 if Z > Z^ (see Fig. [7]). This means that indiscernibility is quantified in the 
discernible domain whereas it is complete in the indiscernible range. This interpretation is 
corroborated by calorimetric measurements. For Z < Z^, we can label 1 or 2, respectively, 
the components of the split wavefunction (@J and ijj n {xj\ — Zd/2) — '0n( a7 j2 + fo/2) 7^ i>n( x j2 — 
ln/2) - ip n {xji + /d/2). Otherwise, for Z > T di , - W 2 ) - ^(^2 + W 2 ) = ^(^2 - 

Zd/2) — ip n ( x ji +/d/2). Consequently, the total molar-entropy of the transition is Z£ In 2 and 
the specific heat jump is AC = 3R/2, as effectively measured.— 

The indiscernibility parameter is also relevant for KsD^^tL^SO^. Dielectric measure- 
ments realize Bloch-states at the center of the Brillouin-zone, via insignificant energy or 
momentum transfer, and variables of statistical-mechanical interest can be calculated from 
only one realization of the thermal pure quantum state.— Protons and deuterons are not 
separate, so the indiscernibility parameter at K is XVi(l p ), where l p = pin + (1 — p)lr>- 
The DIB corresponds to a quantum phase transition at K for p di = 0.74 ± 0.02. The 
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FIG. 7. Temperature laws for the indiscernibility parameter driving the phase transition of 
K3D( 1 _ p )H p (S04)2, according to eqs (jHJ) and (fTUj). Dash: the dividing line between discernibil- 
ity (below) and indiscernibility (above) domains. 



temperature law analogous to (jSj) 



XV»{T < T* 
XVP{T > t p 

XV (T) = 1 



n>l V < p < p di 



} 



Pdi < p < 1 y 



(10) 



accounts for the thermally induced phase transition at for p < p di and for the lack of 
transition for p > Pdi- Here, we assume that and realized through energy transfer 
are separable and independent of p. Visual examination of Fig. [7] shows that for T < Tt, 
XV P is shifted upward by XVi(l p ) —XVi(Iq) (Iq = In) and can be written as: 



1 di ~ 



rpp 

1 di 



1/2-XV^ 

1/2 -xvm Tdl] °^p < p^ 

0; pdi < p < 1 



(11) 



The quasi-linear variation of for p < p^, mirroring that of XVi(l p ), is in reasonably 
good agreement with measurements (see Fig. [8]). In contrast, the Ising model advocated by 
Moritomo et al. yields a curved function with a very steep slope as T c (p) — > (see Fig. 3 
in Ref. 15]), so the extrapolated critical concentration (p c = 0.66 ± 0.04) 15 is smaller than 
Pdi- In practice, dielectric data suggest that it is technically difficult to measure T d i <~ 10 
K with precision, so p c is logically a lower bound for p di . 
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FIG. 8. Solid line: Temperature of the discernible-indiscernible transition for K^D^.^H^SO^, 



according to ec 
Fig. 3 in Ref. 



jlip . •: experimental data (dielectric measurements) estimated from digitized 



la] . DIB: discernible-indiscernible boundary. l p = (1 — p)l^ + pin is the effective 



separation length for the components of the wavefunction. 



According to Fig. [7J the specific heat jump AC(p) for p < pd% is proportional to 



d i 

J TP"(T)dT «T^[l/2-XPi(g], 



(12) 



so 



AC(p) „ SR T^-IV^)] ^ Q< 



(13) 



AC(p) = 0; p di < p < 1 

This is in reasonably good agreement with observations (see Fig. |9]). The shape of AR(p) 
mirrors that of TV°(T) in Fig. [71 It is nearly linear for p 0.4 and the slope goes to zero 
as TP\{l p ) — > lVi(ldi), so there is no discontinuity at the DIB. 



V. Rb 3 H(S0 4 ) 2 



The INS spectrum of the isomorphous crystal of Rb 3 H(S04)2 (see Fig. |2]and Table |V1I|) is 



also consistent with a square- well, apart from the transition observed at 44 cm -1 , instead 
of 57 cm -1 . However, the high background casts doubt on the actual frequency that should 
be reassessed with a more appropriate instrument.— In any case, an OH band in such a 
low frequency range is in conflict with the interpretation of NCS data proposed by Homouz 
et al.— for a single-crystal of Rb3H(S04)2 at 10 K. These authors observed that the mean 
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l p (A) 

0.24 0.22 0.20 0.18 



0.16 



o 



1.5 
1.2 
0.9 
0.6 
0.3 



0.0 




Discernibility s 

(A2/a) 



Indiscernibility 

(A2/a) 



DIB 



0.0 0.2 0.4 0.6 0.8 1.0 
P 



FIG. 9. Solid: specific heat jump for K 3 D( 1 _ p )H /9 (S04)2, according to eq. (TT5j) . Dot dash: the 
ising model, as shown in Ref. [15]. •: experimental data estimated from digitized Fig. 4 in Ref. 
1|]]. DIB: discernible-indiscernible boundary. l p = (1 — p)fo + P^H is the effective separation length 

for the components of the wavefunction. 

proton-momentum values along the crystal axes are rather similar (c^ = 4.60 A -1 , a a = 3.73 
A -1 , <r c * = 4.35 A^ 1 ), from what they inferred that the three OH modes should be in the 
1000-1250 cm -1 range. This is definitively a logical fallacy because the OH modes are not 
parallel to the crystal axes. In fact, we can estimate the mean- momentum for each OH 
mode from the known frequencies. With the square- well model, (pi)] 12 = h/2a « 1 A" 1 .-^ 
For the bending modes observed at ~ 1250 cm -1 and « 1600 cm -1 for both K3H(S04)2 and 
Rb 3 H(S04)2, the harmonic approximation gives {p 2 y )l /2 « 4.39 A" 1 and (p 2 z )l /2 « 5.00 A" 1 . 
Then, the projections along the crystal axes are: 



(p 2 >} /2 cos30° + (p 2 }^m30 c 
(pi) 1 / 2 sin30° + (^)V 2 cos30 c 



3.37 A 1 



«< 1/2 



4.83 A 
4.39 A -1 



> . 



(14) 



(Pl) 1/2 

(Pl) 1/2 

These projections are effectively rather similar although the distribution is markedly 
anisotropic. Furthermore, the low vOH frequency imposes o\, < a a , so we suspect a a 
and o"b were swaped by Homouz et al. Then, taking into account this permutation, nu- 
merics ( I14p are in reasonably good agreement with NCS values: a a — (pD^ 2 = 0.36 A -1 ; 
o h - (p 2 a )l /2 = -0.23 A -1 ; a x - (p 2 c *)\ /2 = -0.04 A -1 . The conflict of interpretation put for- 
ward by Homouz et al. is pointless and there is every reason to suppose that the square-well 
holds for Rb 3 H(S04)2. Our interpretation should be relevant for the dielectric anomaly of 
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Rb 3 D(S0 4 ) 2 observed at 82 K. 



VI. CONCLUSION 

The space group symmetry of K 3 H(S04)2 or K 3 D(S04)2 determined with neutron diffrac- 
tion is A2/a at cryogenic or room temperature and there is no resolvable site splitting for H or 
D. Statistical disorder and antiferroelectricity are ruled out. The anomalies of the dielectric 
constants or heat capacity do not correspond to any structural symmetry breaking. 

At low temperatures, the probability densities for H or D and the INS spectrum of 
K3H(S04)2 are consistent with split-waves representing rigid entities Hli/2— Hri/2 (^h ~ 0.16 
A) or D L1 / 2 — Dri/2 (It> ~ 0.25 A) whose center of mass experiences a square-well. We 
leave open the question as to whether computational chemistry could confirm or not this 
counterintuitive well, as opposed to double-wells for longer bonds or funnel-shaped single- 
wells for shorter bonds. 

Our explanation of the dielectric and calorimetric anomalies is based upon two postulates. 
The first one is that these anomalies correspond to the crossing of a definite boundary for the 
discernibility-indiscernibility degree. To the best of our knowledge, this is unprecedented 
and we ignore whether or not this boundary is deeply rooted in the ground of quantum 
mechanics. The indiscernibility parameter can be thought of as a quantum counterpart 
to an order parameter diverging at the phase transition in classical physics. The second 
postulate is the correlation of the indiscernibility parameter and the dielectric constant, via 
the S— O bonds. The physical ground of this correlation is largely unknown. Our theory is 
scaled by three measured parameters, a, In, Id- It accounts precisely for every observation, 
within the framework of the quantum theory of measurements. The crystal is a macroscopic- 
scale quantum system for which the existence of a thermal pure quantum state is allowed 
by the adiabatic separation of H or D and heavy nuclei, as commonly encountered in similar 
O- • -O bonds. 
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TABLE I. Neutron single crystal diffraction data and structure refinement for K3H(S04)2 and 
K 3 D(S0 4 ) 2 - A = 0.8305 A. Space groups monoclinic Alja with Z = 8. The criterion for used 
reflections was I > 3<r(I). The variance for the last digit is given in parentheses. Refinement on F. 





lS-3rH&U4 / )2 






T<"„TVQn ^„ 
J\3-U^SU4j2 




5 K 


12 K 


293 K 


14 K 


293 K 


a (A) 


9.704(1) 


9.700(1) 


9.806(1) 


9.777(1) 


9.778(1) 


b(A) 


5.639(1) 


5.634(1) 


5.687(1) 


5.674(1) 


5.681(1) 


c(A) 


14.526(1) 


14.524(1) 


14.702(1) 


14.667(1) 


14.701(1) 


P ( ) 


i no nt; /1 \ 

102.95(1) 


1 no to 1 1 \ 

102.73(1) 


1 no no f 1 \ 

102.93(1 ) 


1 no nTi 1 \ 

102.97(1 ) 


1 no it \ 

103.05(1 ) 


v (A 3 ) 


774.5(2) 


774.3(2) 


799.1(2) 


792.9(2) 


795.5(2) 


Measured reflections 


2100 


2385 


2687 


1822 


794 


Independent reflections 


1749 


1373 


1373 


1253 


432 


Used reflections 


1524 


1183 


1781 


999 


404 


Coefficients 


61 


61 


61 


61 


61 


Rint 


0.0002 


0.0015 


0.0002 


0.0004 


0.0006 


R-factor 


0.027 


0.045 


0.024 


0.047 


0.018 


Weighted R-factor 


0.023 


0.051 


0.014 


0.033 


0.020 


Goodness of fit 


1.101 


0.885 


1.077 


1.066 


1.064 


Extinction coefficient 


13.8(5) 


14.8(10) 


9.6(4) 


5.0(3) 


40(1) 
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TABLE II. Atomic positions and isotropic temperature factors for K3H(S04)2 at 5 K (first lines), 
12 K (second lines) and 293 K (third lines). 



Atom 


x/a 


y/b 


z/c 


U(iso)( 


K(l) 


0.2500 


0.75077(11) 


0.0000 


0.0041 




0.2500 


0.7508(2) 


0.0000 


0.0048 




0.2500 


0.75551(13) 


0.0000 


0.0198 


K(2) 


0.34809(5) 


0.73247(8) 


0.30601(3) 


0.0044 




0.34802(9) 


0.73254(16) 


0.30583(7) 


0.0047 




0.34786(5) 


0.73208(8) 


0.30578(4) 


0.0226 


S(l) 


0.03807(6) 


0.2726(1) 


-0.11626(4) 


0.0025 




0.03802(11) 


0.2725(2) 


-0.11641(8) 


0.0024 




0.03807(5) 


0.27302(8) 


-0.11599(3) 


0.0131 


0(1) 


0.09896(3) 


0.47808(5) 


-0.15424(2) 


0.0050 




0.09897(6) 


0.4780(1) 


-0.15423(5) 


0.0055 




0.09819(3) 


0.47694(5) 


-0.15331(2) 


0.0205 


0(2) 


0.10868(3) 


0.05211(5) 


-0.13025(2) 


0.0051 




0.10875(5) 


0.0521(1) 


-0.13021(4) 


0.0053 




0.10805(4) 


0.05494(5) 


-0.13010(2) 


0.0224 


0(3) 


-0.11553(3) 


0.25992(5) 


-0.15300(2) 


0.0050 




-0.11559(6) 


0.2599(1) 


-0.15298(4) 


0.0052 




-O.H38l(3) 


0.26055(6) 


-0.15244(2) 


0.0205 


0(4) 


0.06148(3) 


0.31314(5) 


-0.009127(19) 


0.0052 




0.06141(6) 


0.31321(11) 


-0.00912(4) 


0.0054 




0.06099(4) 


0.31245(7) 


-0.010210(19) 


0.0241 


H(l) 


0.0000 


0.5000 


0.0000 


0.0213 




0.0000 


0.5000 


0.0000 


0.0222 




0.0000 


0.5000 


0.0000 


0.0477 
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TABLE III. Atomic positions and isotropic temperature factors for K3D(S04)2 at 14 K (first lines) 
293 K (second lines). 



Atom 


x/a 


y/b 


z/c 


U(iso)( 


K(l) 


0.2500 


0.7517(3) 


0.0000 


0.0042 




0.2500 


0.7559(2) 


0.0000 


0.0217 


K(2) 


0.34821(12) 


0.7343(2) 


0.30569(9) 


0.0052 




0.34804(11) 


0.73270(17) 


0.30574(8) 


0.0242 


S(l) 


0.03846(15) 


0.2704(3) 


-0.1164(1) 


0.0027 




0.03823(13) 


0.2727(2) 


-0.11604(8) 


0.0142 


0(1) 


0.09942(8) 


0.47698(14) 


-0.15382(5) 


0.0048 




0.09827(6) 


0.47659(11) 


-0.15307(4) 


0.0219 


0(2) 


0.10852(9) 


0.05043(15) 


-0.13138(5) 


0.0057 




0.10803(6) 


0.05459(12) 


-0.13039(5) 


0.0238 


0(3) 


-0.11496(8) 


0.25884(15) 


-0.15246(5) 


0.0050 




-0.11348(7) 


0.26004(11) 


-0.15222(4) 


0.0221 


0(4) 


0.06275(8) 


0.30902(15) 


-0.00884(5) 


0.0055 




0.06161(7) 


0.31036(16) 


-0.01022(4) 


0.0261 


D(l) 


0.0000 


0.5000 


0.0000 


0.0396 




0.0000 


0.5000 


0.0000 


0.0533 
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TABLE IV. Thermal parameters in A 2 units for K3H(S04)2 at 5 K (first lines), 12 K (second 
lines) and 293 K (third lines). The thermal parameters account for the variation of the contri- 
bution of each atom to Bragg's peak intensities through the thermal factor T at depending on 
the reciprocal lattice parameters a*, b* , c*, and unit cell indexes in reciprocal space h,k,l, as 
T at = exp[-27r 2 ([/ 1 a */ i 2 a* 2 + U%k 2 b* 2 + Ugl 2 c* 2 + 2U$ha*kb* + 2Ugkb*lc* + U${lc*ha*)\. 



Atom 


Un 


u 22 


u 33 


u 23 


U13 


U12 


K(l) 


0.0043(2) 


0.0048(2) 


0.0031(2) 


0.0000 


0.00068(16) 


0.0000 




0.0037(5) 


0.0046(5) 


0.0062(6) 


0.0000 


0.0010(4) 


0.0000 




0.0184(2) 


0.0242(2) 


0.0165(2) 


0.0000 


0.00314(16) 


0.0000 


K(2) 


0.00446(16) 


0.00402(16) 


0.00458(15) 


-0.00050(12) 


0.00104(12) 


-0.00024(12) 




0.0042(4) 


0.0037(4) 


0.0063(4) 


-0.0008(3) 


0.0010(3) 


-0.0004(2) 




0.01899(15) 


0.02015(16) 


0.02811(19) 


-0.00470(15) 


0.00395(14) 


-0.00035(14) 


S(l) 


0.00272(19) 


0.00217(18) 


0.00252(18) 


0.00001(14) 


0.00057(14) 


0.00014(15) 




0.0031(4) 


0.0023(5) 


0.0021(5) 


-0.0004(3) 


0.0009(3) 


0.0001(3) 




0.01257(15) 


0.01287(15) 


0.01369(16) 


0.00044(12) 


0.00284(12) 


0.00021(13) 


0(1) 


0.00599(12) 


0.00410(11) 


0.00525(11) 


0.00060(8) 


0.00187(8) 


-0.00108(8) 




0.0058(2) 


0.0036(3) 


0.0073(3) 


0.00072(17) 


0.00221(18) 


-0.00122(15) 




0.02433(13) 


0.0168(1) 


0.02191(11) 


0.00121(8) 


0.00822(9) 


-0.00460(8) 


0(2) 


0.00578(11) 


0.00396(11) 


0.00583(11) 


-0.00002(8) 


0.00169(8) 


0.00145(8) 




0.0054(2) 


0.0033(3) 


0.0074(3) 


0.00026(17) 


0.00179(19) 


0.00153(15) 




0.02339(13) 


0.01679(11) 


0.02721(13) 


0.00138(9) 


0.00609(9) 


0.00638(9) 


0(3) 


0.00338(11) 


0.00607(11) 


0.0053(1) 


-0.00023(8) 


0.00057(7) 


-0.00014(8) 




0.0034(2) 


0.0056(3) 


0.0063(3) 


-0.00033(17) 


0.00058(17) 


-0.00018(14) 




0.01303(9) 


0.02581(12) 


0.02203(11) 


-0.0023(1) 


0.00271(7) 


-0.00040(9) 


0(4) 


0.00666(12) 


0.00599(12) 


0.0029(1) 


-0.00011(7) 


0.00090(8) 


0.00021(8) 




0.0065(2) 


0.0056(3) 


0.0040(3) 


-0.00022(17) 


0.00078(17) 


0.00022(15) 




0.02840(14) 


0.03032(14) 


0.0127(1) 


0.00036(9) 


0.00304(9) 


0.0003(1) 


H(l) 


0.0209 


0.0304 


0.0129 


-0.0019 


0.0041 


-0.0082 




0.0204 


0.0282 


0.0178 


-0.0017 


0.0041 


-0.0073 




0.0457 


0.0765 


0.0212 


-0.0045 


0.0081 


-0.0282 
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TABLE V. Thermal parameters in A 2 units for K3D(S04)2 at 14 K (first lines) and 293 K (second 
lines) . The thermal parameters account for the variation of the contribution of each atom to Bragg's 
peak intensities through the thermal factor T at depending on the reciprocal lattice parameters a*, 
b*, c*, and unit cell indexes in reciprocal space h,k,l, as T at = exp[-27r 2 (C/f 1 */i 2 a* 2 + Ugffib* 2 + 
Ugl 2 c* 2 + 2Uft,ha*kb* + 2Ugkb*lc* + U$lc*ha*)]. 



Atom 


Un 


u 22 


u 33 


u 23 


U13 


U12 


K(l) 


0.0037(6) 


0.0056(7) 


0.0032(6) 


0.0000 


0.0006(5) 


0.0000 




0.0170(9) 


0.028(1) 


0.0196(8) 


0.0000 


0.0030(7) 


0.0000 


K(2) 


0.0043(4) 


0.0041(5) 


0.0071(4) 


-0.0015(4) 


0.0012(3) 


-0.0003(4) 




0.0174(7) 


0.0240(7) 


0.0301(7) 


-0.0057(5) 


0.0033(5) 


-0.0006(4) 


S(l) 


0.0023(5) 


0.0027(6) 


0.0026(5) 


0.0008(5) 


-0.0002(4) 


-0.0004(5) 




0.0109(7) 


0.0155(7) 


0.0161(7) 


0.0019(5) 


0.0026(5) 


0.0004(4) 


0(1) 


0.0057(3) 


0.0030(3) 


0.0060(3) 


0.0006(2) 


0.0020(2) 


-0.0012(2) 




0.0227(4) 


0.0200(5) 


0.0243(5) 


0.0010(3) 


0.0079(3) 


-0.0051(3) 


0(2) 


0.0063(3) 


0.0040(3) 


0.0070(3) 


0.0002(3) 


0.0020(2) 


0.0016(3) 




0.0216(4) 


0.0205(5) 


0.0293(4) 


0.0011(3) 


0.0059(3) 


0.0072(3) 


0(3) 


0.0033(3) 


0.0058(3) 


0.0059(3) 


-0.0012(3) 


0.0009(2) 


-0.0004(3) 




0.0115(4) 


0.0298(5) 


0.0243(5) 


-0.0028(3) 


0.0029(3) 


-0.0009(2) 


0(4) 


0.0075(3) 


0.0055(3) 


0.0036(3) 


-0.0013(3) 


0.0016(2) 


0.0018(3) 




0.0280(5) 


0.0335(5) 


0.0158(5) 


0.0003(3) 


0.0028(3) 


0.0003(4) 


D(l) 


0.035(1) 


0.0728(16) 


0.0112(6) 


-0.0053(7) 


0.0066(5) 


-0.035(1) 




0.0464(7) 


0.0907(11) 


0.0222(6) 


-0.0052(6) 


0.0066(5) 


-0.0370(8) 
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TABLE VI. Bond lengths in A units of the SO4 entities and their expansion A with temperature. 







K 3 H(S0 4 ) 2 






K 3 D(S0 4 ) 2 










5 K 


293 K 




14 K 


293 K 


AR» Q 


Sl- 


-01 


1.4636(6) 


1.4623(5) 


0.0013(11) 


1.4754(17) 


1.4585(14) 


0.0169(31) 


Sl- 


-02 


1.4555(6) 


1.4545(6) 


0.0010(12) 


1.4637(17) 


1.4527(13) 


0.0110(30) 


Sl- 


-03 


1.4680(6) 


1.4674(5) 


0.0006(11) 


1.4758(16) 


1.4602(15) 


0.0156(31) 


Sl- 


-04 


1.5385(6) 


1.5372(5) 


0.0013(11) 


1.5577(16) 


1.5354(14) 


0.0223(30) 



TABLE VII. Diagonal thermal factors in A 2 units: U H / Dx (T) parallel to O- • O and U H / Dz (T) 
perpendicular to O- • O and parallel to (a, b). 



U Hx (5) = 0.0351 


U Ux (12) 


= 0.0326 


U Hx (293) 


= 0.0932 




U Dx (U) 


= 0.0937 




= 0.1117 


U Hz (5) = 0.0162 


Ch*(12) 


= 0.0160 


f/ H2 (293) 


= 0.0290 




U Dz (U) 


= 0.0141 


U Dz (293) 


= 0.0254 



TABLE VIII. Computed half- widths at half-height in A units of the probability densities: A R / Dx (T) 
parallel to O- • O; A H / Dz (T) perpendicular to O- • O and parallel to (a, 6); A n / Dy (T) parallel to 
c*. 



A Hx (5) = 0.156 


A Hz (5) = 0.106 


A Hy (5) = 0.095 


A Hx (12) = 0.150 


A Hz (12) = 0.105 


A Hy (12) =0.111 


A Hx (293) = 0.254 


A Hz (293) = 0.145 


A Hy (293) = 0.142 


A Dx (14) = 0.255 


A Dz (14) = 0.099 


A Dy (14) = 0.088 


A Dx (293) = 0.278 


A Dz (293) = 0.133 


A Dy (293) = 0.124 



24 



TABLE IX. Tentative assignment scheme of the OH stretching transitions observed with inelastic 
neutron scattering at 20 K, from Ref. [lj], and comparison with the calculated transitions for a 
bare proton in an infinite square well whose width is a = 0.47 A 







1 - 2 


1 - 3 


1-4 


1 - 5 


1 - 6 


K 3 H(S0 4 ) 2 


Obs. (cm -1 ) 


57 ±3 


145 ±5 


280 ± 10 


445 ± 15 


680 ± 20 




Calc. (cm -1 ) 


57 ±3 


152 ±8 


285 ± 15 


456 ± 24 


665 ± 35 


Rb 3 H(S0 4 ) 2 


Obs. (cm -1 ) 


44±? 


145 ±5 


260 ± 10 


440 ± 15 


680 ± 20 
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Projections of dimer entities onto (a, b) planes. Thermal ellipsoids represent 
50% of the probability density for nuclei. Sulfur is yellow, oxygen is red, 
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INS spectra of K 3 H(S04)2 and Rb3H(SOzi)2 crystal powders at 20 K, after 



ref. 14J. The bands at 600 cm -1 were assigned to SO4 entities 

3 Eigenstates of the square-well along the OH stretching coordinate 

4 Solid lines: Computed probability densities along the OH stretching coordi- 
nate at 5 K (bottom) and 293 K (top) for a rigid H1/2— H1/2 entity. The length 
is /h = 0.16 A and the center of mass experiences a square- well whose width 
is 0.47 A. Dashed line: probability density for each half proton. +: Gaussian 
profile whose variance is U^(T)/2 (see Table VII) 

5 Solid lines: Computed probability densities along the OD stretching coordi- 
nate at 14 K, 85 K, and 293 K for a rigid T>li/2— ~Dri/2 entity. The length is 
Id = 0.25 A and the center of mass experiences a box whose width is 0.47 A. 
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6 (a) The indiscernibility degree of the ground state, X^Z), as a function of 
the separation of the split-wave. Zh = 0.16 A, T\ = 0.60 for K 3 H(S0 4 ) 2 . 
Z D = 0.25 A, Xi = 0.27 for K 3 D(S0 4 ) 2 . hi = 0.183 A, l ldi = 0.50 at the 
discernible- indiscernible boundary, (b) The transition parameter XV = I± if 
I < ki, TP = \\il>l dl 

7 Temperature laws for the indiscernibility parameter driving the phase transi- 
tion of K 3 D(i_p)Hp(S04)2, according to eqs (8) and (10). Dash: the dividing 
line between discernibility (below) and indiscernibility (above) domains 
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Solid line: Temperature of the discernible-indiscernible transition for K 3 D( 1 _ /3 )H P (S04)2, 
according to eq. (11). •: experimental data (dielectric measurements) es- 
timated from digitized Fig. 3 in Ref. 15) . DIB: discernible- indiscernible 
boundary. l p = (1 — p)l^ + pZn is the effective separation length for the 

components of the wavefunction 14 

Solid: specific heat jump for K 3 D(i_p)H p (S04)2, according to eq. (13). Dot 



dash: the Ising model, as shown in Ref. 



15) . •: experimental data estimated 



from digitized Fig. 4 in Ref. 151 ] . DIB: discernible-indiscernible boundary. 



l p = (1 — p)Z D + p/ H is the effective separation length for the components of 



the wavefunction 15 
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